Behavior of test particles in the plasma sheath upon local bias and plasma switching 
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Equilibrating gravitation by electric forces, microparticles can be confined in the plasma sheath 
above suitably biased local electrodes. Their position depends on the detailed structure of the 
plasma sheath and on the charge that the particles acquire in the surrounding plasma, that is by 
the electron and ion currents towards it. Bias switching experiments reveal how the charge and 
equilibrium position of the microparticle change upon altered sheath conditions. We observe similar 
particle behaviors also in the afterglow of the discharge for a persisting bias voltage on the electrode: 
damped oscillation into a new equilibrium or (accelerated) fall according to the bias. 

PACS numbers: 52.27.Lw, 52. 80. Pi, 52.65.Rr, 52.40.Kh 



I. INTRODUCTION 

At the powered electrode of an asymmetric rf- 
discharge, self biasing gives rise to a large potential drop, 
wide plasma sheath and strong electric fields. In contrast, 
the plasma sheath in front of grounded or additionally bi- 
ased electrodes is less pronounced and thinner. The prop- 
erties of such sheaths are, nevertheless, of fundamental 
interest in view of plasma processing of surfaces^ and 
plasma chemistryP Selective biasing of the substrate can 
be used to optimize the electron and ion impact energies 
and directions.^ 

Monitoring plasma characteristics in the sheaths 
is challenging, since common diagnostics, such as 
Langmuir-probes usually fail. This is because the macro- 
scopic tip of a Langmuir-probe significantly alters the 
local plasma characteristics and cannot resolve density 
and potential gradients in the narrow sheath with ade- 
quate precision. In order to overcome these limitations, 
alternative diagnostics like optical spectroscopy^ or mi- 
croparticles as probes^] are used. Those provide an effi- 
cient, nearly non-invasive tool for the local diagnostic of 
the plasma parameters. In view of their practical use, a 
calibration is necessary to discern their own charge dy- 
namics from their capability to reflect the properties of 
the surrounding plasma. 

Upon immersion in a plasma, microparticles get 
charged by electrons and ions that accumulate on their 
surface. Since the electron mobility exceeds those of the 
ions, the particle charge is in general negative, depending 
on the local conditions of the surrounding plasma. Al- 
beit, it has been stressed that under certain conditions, 
such as in a discharge afterglow, microparticles in the 
bulk plasma can even acquire positive charges. 9 While 
the plasma is operational, particle recharging takes place 
on longer time scale than the almost instantaneous re- 
action of the plasma to slightly modified external con- 
ditions. Then, the particle charge constantly remains 
in equilibrium with the surrounding electrons and ions. 



Drastic changes in the operating conditions, however, 
may perturb the time scale ordering and hinder an equi- 
libration of the particle charge. 

In this work, we investigate microparticles which 
are captured in the plasma sheath and their reaction 
to changes in the operating conditions of the plasma. 
Thereby, we focus on switching off the plasma while the 
lower discharge electrode is grounded or additionally bi- 
ased as well as changing the bias voltage during operation 
of the plasma. The experimental results are discussed in 
view of a theoretical model that combines basic equations 
with input data from particle-in-cell (PIC) simulations. 



II. EXPERIMENTAL SETUP 

The PULVA-INP setup 10 consists of an asymmetric, 
capacitively coupled rf-plasma in argon, working at neu- 
tral gas pressures pAr from 0.1 Pa to 100 Pa. The rf- 
power (P r f [W] = 5 — 100) is supplied by the upper, pow- 
ered electrode at a frequency of v x { = 13.56 MHz and 
amplitudes $ r f up to 1000 V. In dependence on the ex- 
ternal parameters P r f and £>An the obtained character- 
istics for the pristine argon plasma are electron densi- 
ties n e [cm -3 ] = 10 9 — 10 11 , electron energies k-BT e [eV] — 
0.8—2.8 and plasma potentials with respect to the ground 
of 3> P [V] = 20 — 30. 11 The overall plasma characteristics 
are monitored by Langmuir probe and plasma monitor 
measurements. 

The key feature of the experimental setup is the lower, 
so called 'adaptive' electrode (AE). It consists of 101 
square electrode segments (pixels) with a linear exten- 
sion of 6.6 mm each, separated by thin (0.4 mm) isolating 
gaps (see Fig. [I]). In addition, four larger segments fit the 
pixel geometry to the surrounding ring and ground shield. 
All 105 electrode pixels can be biased individually or in 
groups by an external dc- voltage of up to $bias = ±100 V. 
The selective application of bias voltages to some pix- 
els allows for studying spatial and temporal changes of 



2 




□□□■aEL 

□□□r ~ 

□nnn 
□□□ 
□□□ 




FIG. 1: (Color online) Setup of the adaptive electrode (only 
one quarter is shown). From the center outward: square 
pixels, fitting segments (green), ring (light grey) and ground 
shield (dark grey). The spacings between all constituents are 
isolating. 



the plasma sheathP^ The small extension of the pixels 
as compared to the remaining grounded electrode guar- 
antees that the applied bias only locally influences the 
plasma sheath but leaves the overall discharge conditions 
unaltered. 

All experiments in this work refer to a plasma op- 
erating power P r f = 10W and neutral gas pressure 
PAr = 5 Pa unless stated differently. The microparti- 
cles used as probes are melamine-formaldehyde (MF) 
spheres with diameter d = 9.6 jam. With a mass den- 
sity of p m = 1.51 x 10 3 kgm -3 for MF this results in a 
particle mass of m = 7 x 10 -13 kg. Illuminated by a laser 
fan (532 nm), their positions are recorded by a fast CCD 
camera (up to 2000 frames per second). 



III. RESULTS AND DISCUSSION 

In this work we report on three different experiments. 
Each one describes a particular way in which the equi- 
librium conditions for a levitating particle in the plasma 
sheath are perturbed. The initial position of the particle 
is determined by an equilibrium between gravitational 
and electrostatic forces. Neutral drag and thermophore- 
sis effects are negligible. The effect of ion drag forces will 
be discussed in the experiments for which it is relevant. 



A. Plasma afterglow — unbiased AE 

In this experiment we start from an initial configu- 
ration where the MF-particle is trapped above the un- 
biased center pixel by a confining potential of — 5V on 
the surrounding pixels of the AE. At pAr = 5 Pa, the 
sheath width is about 3 mm and the equilibrium posi- 
tion zq = 1.8 mm. Switching off the plasma, the electro- 




FIG. 2: (Color online) Main panel: Measured vertical particle 
position after switching off the plasma for different neutral gas 
pressure pAr (boxes). Theoretical trajectories for freely falling 
[retarded by friction according to |l])] particles are given by 
dashed-dotted [solid] lines. Inset: Crosses give the friction co- 
efficient P/m as extracted from fitting |l]) to the experimental 
data. The dotted line is a guide to the eye. 



static force no longer compensates gravity and the par- 
ticle drops onto the AE. As compared to a freely falling 
particle, the drop is markedly retarded (see Fig. [2|. 

Along with the residual electric field, which decays over 
some /is, also the directed ion current towards the AE 
vanishes and ion drag can be neglected. The only forces 
relevant on the time scale of the fall (ms) are gravitation 
F g = —mge z and neutral drag force F n = — /3z = —f3ze z . 
Here the damping constant f3 depends on the local den- 
sity of the background gas, which we assume to be con- 
stant in time. Then the equation of motion for the par- 
ticle reads z + ^z + g = 0. The solution of this equa- 
tion, subject to the initial conditions z(t = 0) = and 
z(t = 0) = zq is given by 



Z(r) Zq 



m 2 g 

P 2 



(1 



(1) 



where we introduced the normalized time r = fit/m. De- 
pending on the argon gas pressure, the values of f3/m 
obtained by a least squares fit of to the experimental 
data are in the range of 20s _1 to50s _1 . This so called 
the Epstein friction coefficient!^ can be related to the 
neutral gas pressure by 



/3 
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where is the thermal velocity of the argon gas atoms 
and S is a parameter in the order of 1 — 1.5 account- 
ing for how the gas atoms are deflected from the particle 
surface!^ Assuming a diffuse reflection, we obtain a the- 
oretical value of /3/m ~ 30 s -1 for pa t = 20 Pa at room 
temperature, which is in good accordance with the values 



3 







Udc = 


-100V 




Udc = 


-80V 




U d c = 


-60V 




Udc- 


-40V 




U d c = 


-5V 



60 






U d c = 


-80V 




Udc- 


-40V 




U d c = 


-20V 




Udc- 


OV 




U d c = 


+5V 



40 

t[ms] 



FIG. 3: (Color online) Measured vertical particle position 
after switching off the plasma while retaining dc-biases on se- 
lected segments of the AE. Left panel: Bias voltages Uae = 
(Ojt/dcjO), where the three numbers refer to the pixel bias- 
ing from the center of the AE outward - (light blue, vio- 
let, dark blue) with pixel colors referring to Fig. [I] The 
rest of the AE is grounded. Right panel: Bias voltages 
U ae = (C7dc,-100V,-100V). Time resolution of the mea- 
surement is 0.5 ms and symbols are given to distinguish the 
curves. Error bars for the positions are comparable to the 
symbol size. 



determined by the measurement. Furthermore, to de- 
scribe our experimental data it is not necessary to adapt 
/3 to the decaying neutral gas density in the discharge 
afterglow, as discussed in Ref. 9. Within our accuracy a 
constant value of j3 is sufficient. 

Due to the vanishing electric field this experiment does 
not provide any information on the value or variation of 
the particle charge. In order to obtain information on 
this aspect, we have to resort to a slightly modified setup, 
retaining a static dc-bias at some pixels of the AE. 



B. Plasma afterglow — biased AE 

Extending the previous experiment, we start from var- 
ious bias configurations of the AE by which we initially 
trap a particle. Clearly, the biased pixels locally warp 
the plasma sheath. Thus, the equilibrium positions (and 
presumably also the charges) of the particles differ for 
each case. The configurations for which the results are 
given in Fig. [3] differ in the number of biased pixels - 
(3 x 3) for the left and (5 x 5) for the right panel, respec- 
tively. The applied bias voltages exhaust the available 
range provided by the AE, resulting in weak to strong 
confinement. For weak confinement {\U^c\ < 60 V) we 
recover in the left panel of Fig. [3] the falling behavior 
observed in the previous experiment. Only for stronger 
negative bias (|?7dc| > 60 V) the persisting electric field 
is strong enough to keep the particle hovering above the 
AE. Note that the lengths of the recorded time inter- 
vals are limited by the particle horizontally leaving the 
camera focus. Their levitating state persists over several 
seconds. The particle trajectories in the right panel of 



20 L U AE [V] = (0, -100, 0) j 

! i J 




60 
50 
40 
30 
20 

\ 10 



-20 -10 



10 20 



O 
> 

i_lT 




20[ U AE [V] = (-80,-100,-100) 

E 

E.10I 



-20 -10 10 20 

x [mm] 



II : — 



40 F 

20 

E 
o 

> -20 
LU N -40 
-60 



plasma on 
plasma off 



5 10 15 20 25 
z [mm] 





■\-v 



— ■ U AE [V] = ( 0, -5, 0) 

— U AE [V] = ( 0, -100, 0) 

■ - U AE [V] = ( 0,-100,-100) 

— U AE [V] = (-80,-100,-100) 



15 

z [mm] 



20 



25 



FIG. 4: (Color online) Left panel: Vertical cut through the 
electrostatic potential distribution above the AE in absence 
of any plasma. The shown region depicts only a fraction of 
the calculation grid (n x x n z = 1400 x 2000 grid points cor- 
responding to 7 cm x 10 cm, x > 0). Right panel: Vertical 
component of the electric field on the symmetry axis as ex- 
tracted from the potential data. In the inset the electric field 
for C7ae[V] = (0,-100,0) before and after switching off the 
plasma is given. The data for operating plasma result from a 
PIC simulation. The diamond indicates the measured initial 
equilibrium position of the particle. 



Fig. [3] are shifted to markedly larger distances from the 
AE. This reflects the more pronounced widening of the 
plasma sheath by the larger area of strongly biased pix- 
els. As long as the surrounding pixels are on a sufficiently 
negative bias, the particle levitation even persists for a 
moderate positive voltage on the center pixel. For a large 
negative bias (U& c = —80 V), however, switching off the 
plasma causes an abrupt drop of the particle, which is 
even accelerated as compared to a free fall. 

In order to explain the observed behavior, a first idea 
would be to assume that on the observed time scales the 
plasma has completely decayed, i.e. all electrons and ions 
have already recombined. Then, the particle trajectory 
corresponds to the motion of a charged sphere between 
two parallel plates, that are the powered electrode and 
the AE. Using the superposition principle, we may calcu- 
late the electrostatic force acting on the particle in two 
steps. First, the self-force for a charged particle between 
two parallel, grounded planes is calculated analytically 
by means of an infinite series of image charges^ (see 
Appendix A). This contribution is, however, negligible 
for distances z > 0.1mm. Second, we numerically solve 
the boundary value problem and calculate the potential 
due to the biased electrode pixels. To this end, we reduce 
the simulation volume to a two-dimensional, cylindrically 
symmetric domain by smearing out the individual pixels 
to rings, and use routines from FISHPACKP In vertical 
direction the simulation volume covers the whole range 
between both electrodes. Radially, the simulation do- 
main is such large that the boundaries do not affect the 
results near the center of the discharge. At r max we used 
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vanishing radial field boundary conditions. 

Due to rotational symmetry, the biasing of outer ly- 
ing pixels influences the electrostatic potential above the 
center pixel to the point of generating a region with neg- 
ative potential at finite z (left panel of Fig.[4|. Above an 
unbiased center pixel this gives rise to a potential struc- 
ture with a minimum at zq. For strong negative biasing 
of the center pixel, U A e[V] = (-80, -100, -100), the po- 
tential along the symmetry axis monotonically increases 
with z. The sign change of the electric field at zq implies 
that a negatively charged particle may only remain hov- 
ering for z > zq since there gravity and electrostatic force 
are antiparallel (right panel of Fig. [4| . As compared to 
the case of operating plasma, the electric field after the 
switching off is markedly reduced (inset of the right panel 
of Fig. [3]). Therefore, keeping up the force equilibrium 
requires an increasing of the particle charge. Whether 
the particle remains hovering or drops down onto the 
AE depends on its ability to acquire enough negative 
charges on a sufficiently short time scale. For systems 
with markedly higher operating pressures the relevant 
stages of a plasma decay have been identified by Couedel 
et al.^, focussing on the bulk properties. The proper- 
ties of the sheath in the afterglow of a pulsed inductively 
coupled rf-plasma in hydrogen have been addressed by 
Osiac et alP^ For the locally biased plasma sheath the 
following aspect comes into play. After the switching off 
the periodic flooding of the sheath with electrons dur- 
ing each rf-period ceases, but electrons and ions contin- 
uously diffuse from the bulk into the sheath. Here, the 
flow of both species is controlled by the applied local bias 
voltages. The charge equilibrium in the sheath will vary 
only slightly as long as enough electrons and ions from 
the bulk enter the sheath to compensate losses at the 
wall. A particle adapts its charge to the altered electron- 
and ion density in surrounding on time scales of some 
1 — 10 /isP^ Since this time scale is faster than the changes 
of the plasma sheath in the discharge afterglow, the par- 
ticle charge constantly remains in equilibrium with the 
charge densities in its surrounding. In order to explain 
the accelerated fall for U AE [V] = (-80, -100, -100) con- 
sider that due to the strong negative bias no electrons 
at all will be able to penetrate into the sheath above the 
center pixel anymore. Therefore, the exclusive presence 
of ions in the surroundings of the particle will drasti- 
cally reduce its negative charge. One might even have to 
consider a positive recharging since gravity and ion drag 
f orc ^!8jt20] are weak to account for the observed be- 
havior alone. The individual contributions to the total 
force causing the accelerated fall will be analyzed for the 
bias switching experiment in the next section. 



C. Persisting plasma — bias switching 

While switching off the plasma certainly alters the 
equilibrium conditions for the particle drastically, the 
changes in densities and fields occur on a moderate time 
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FIG. 5: (Color online) Measured vertical particle position 
upon bias switching while the plasma is kept operational. The 
initial confining potential is ?7ae[V] = (0,-5, —50), which is 
switched to the final configuration Uae = (Udc, Udc, —50 V). 
In the inset the time scale is magnified to resolve the dynamics 
of those particles that drop onto the AE. For comparison the 
trajectory of a freely falling particle is given. 



scale. This is in contrast to the situation where the AE- 
pixel biasing is switched during operation of the plasma. 
Here, the local plasma sheath adapts within 1 — 5 re- 
cycles to the new biasing configuration and thus has an 
instantaneous impact on the charge and force equilibrium 
of the particle. In principle, a particle may be trapped 
into a stable equilibrium position for all indicated bias- 
ing configurations in Fig. [5j provided that the relaxation 
is performed smoothly. Upon an abrupt change from an 
initial configuration, however, the particle drops onto the 
AE outside a rather narrow region of bias variations (see 
inset of Fig. [5| . This holds for either positive or negative 
voltages. Interestingly, again a strong negative bias leads 
to an acceleration of the particle towards the AE which 
is not explicable solely by gravity. For moderate changes 
in the biasing configuration the particle reaches its new 
equilibrium position through a damped oscillatory mo- 
tion. 

In order to describe the different behaviors quantita- 
tively, we performed a PIC-simulation of the plasma dis- 
charge. Special emphasis is put on the local modifica- 
tions of the plasma sheath by the AE-pixels. Details on 
the simulation can be found in Appendix [Bj 

The initial AE-configuration in Fig. [5] only marginally 
disturbs the plasma sheath above the center pixel. Its 
primary purpose is to keep the particle in the camera fo- 
cus. As can be seen from the simulation results in Fig.[6j 
both the electron and ion densities near the symmetry 
axis agree with those in the unperturbed sheath far away 
from the biased pixels. In this respect, the results for 
^ae[V] = (0,-5,-50) in Fig. 7] are representative for 
an unperturbed sheath. While the two-dimensional color 
plots in Fig. [6] provide a general overview of the sheath 
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FIG. 6: (Color online) PIC-results for the time-averaged po- 
tential as well as electron and ion density. Only part of the 
total simulation volume is shown. The AE-pixels are indi- 
cated by hatched rectangles and the biasing corresponds to 
the initial configuration in Fig. [5] In the upper panel Ri — R3 
indicate the diagnostic regions for which the local electron 
energy probability distributions are given in Fig. [8] 
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FIG. 7: Influence of the pixel biasing on electric field as well 
as electron and ion density on the symmetry axis of the dis- 
charge. Results obtained by PIC simulation. 



structure above the AE-pixels, Fig. [7] allows for a quanti- 
tative analysis on the symmetry axis. The time-averaged 
electron density is markedly reduced as compared to the 
ion density. Electrons may overcome the potential barrier 
of the plasma potential and enter the sheath only during 
a short fraction of an rf-period. This occurs, when the 
potential on the powered electrode is sufficiently negative 
to supply enough kinetic energy to the electrons. Upon 
additional negative biasing of an AE-pixel, the necessary 
energy for crossing the sheath increases. This results 
in a further reduction of the electron density since for 
most electrons the sheath becomes impenetrable. Hav- 
ing lost their kinetic energy already after some part of 
the sheath, they remain trapped in the plasma. Hence, 
above the negatively biased pixel, parts of the sheath 
may be completely deprived of electrons. Because of the 
continuity equation and acceleration of the ions towards 
the AE, their (time- averaged) density also decreases in 
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FIG. 8: (Color online) Left panel: PIC results for the local 
electron energy probability function g p (E) in the bulk and 
the sheath of the discharge. Right panel: Changes in g p (E) 
above the center AE-pixel induced by different biasing. Away 
from the biased pixels, the EEPF in region R3 agrees for all 
AE configurations. 



the sheath. Their larger mass makes them less suscep- 
tible to further acceleration by additional biasing and 
therefore all ion densities nearly coincide. In contrast to 
switching off the plasma, which leads to wider final po- 
tential structures and weaker fields, a negative bias may 
drastically enhance the electric field in the sheath. This 
enhancement results from the dominance of the increased 
potential difference over the effect of sheath widening. 

The density profile does not contain any information 
on the energy of the electrons in the sheath. There- 
fore, we calculate the electron energy probability distri- 
bution g p (E)Wln Fig.Mg p (E) is given for three regions of 
the discharge, representing bulk, unperturbed sheath and 
sheath above a biased electrode pixel. The scales used in 
Fig. [8] are chosen such that a Maxwellian distribution 
results in a straight line. Then the electron tempera- 
ture is proportional to the inverse slope and the crossing 
with the ordinate gives the species density. The obtained 
g p (E) for the bulk closely resembles a bi-Maxwellian dis- 
tribution with a larger fraction of cold and and a smaller 
fraction of higher energetic electrons. The latter one is 
due to those electrons that oscillate between the elec- 
trodes of the discharge. In the sheath, the density is re- 
duced for all energies, but most pronounced in the high 
energy tail. This is due to the retardation caused by the 
potential difference the electrons encounter when cross- 
ing the sheath. As compared to the unperturbed sheath, 
g p (E) above the center pixel clearly reflects the influence 
of the additional AE-biasing (right panel of Fig. [8]). For 
energies up to the plasma potential, the density reduction 
(enhancement) for negative (positive) bias is uniform for 
all energies. The strong fluctuations at larger energies 
are due to poor statistics since only a tiny fraction of 
all electrons have energies in this range. Irrespective of 
the biasing on the center pixel, the confinement poten- 
tial on the third shell of pixels increases the density in the 
high energy tail as compared to the unperturbed sheath 
through a channeling effect. 

Explaining the accelerated fall of the particle quanti- 
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FIG. 9: (Color online) Upper left panel: Total force F^ al 
on the particle as extracted from the experimental trajectory. 
Individual contributions to F^ al are gravitational force F g , 
ion drag force Fi and electric field force Ff. The shown F\ 
for a particle with q — — 10 5 eo is an upper bound estimate of 
this contribution. Ff is calculated for various positive particle 
charges using the PIC results for the electric field. Upper right 
panel: Collection radius b c and Coulomb scattering radius b s 
for different particle charges q as derived from the PIC data 
in the lower panel. Curves and legend correspond from top to 
bottom. Lower panel: PIC results for ion density rii, ion ther- 
mal velocity vl h and ion drift velocity Vi above the center pixel 
of the AE for bias configuration Uae\V\ = (-20, -20, -50). 



tatively requires a careful analysis of all relevant forces. 
In order to extract the total force on the particle from 
the experimental trajectory, we performed a polynomial 
fit to the measured data (inset of Fig. [5| and calculated 
its second derivative. Due to the limited number of data 
points, the curvature data depend slightly on the order 
of the used polynomial, but they agree within 5% in the 
relevant time interval. Combining trajectory and time 
dependent force data, the total force on the particle is 
given as a function of distance above the AE in the up- 
per panel of Fig. [9] Gravity is responsible for only about 
10% of the observed particle acceleration. Also the ion 
drag force, which presumably becomes enhanced by the 
negative pixel biasing, contributes to the acceleration. It 
is given by 



Fi = miniViVi(iib 2 c + 7r6 2 In A) 



(3) 



where mi, ni, v\ are the ion mass, density and drift ve- 
locity. Accounting for the thermal motion of the ions, we 
define the effective velocity vf = vf + (v* 11 ) 2 , where vf 1 
is the thermal ion velocity. In (J3|, the two contributions 
to the ion drag force correspond to momentum transfer 
due to collection of ions (index c) and Coulomb scatter- 
ing (index s). The radii in the associated cross sections 
are given by^ 
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where is the potential of the dust particle, and 
eo the permittivity of vacuum. The divergence of the 
Coulomb scattering cross section is circumvented by us- 
ing a su itable cu toff beyond which passing ions are 
neglecte d 4 * 18 * 20 * 21 -! The shielding effect of the plasma 
suggests 20 to use the electron Debye length Ad = 
[(eo/cBT , e )/(n e eo)] 1 ^ 2 as cutoff, leading to the approximate 
Coulomb logarithm In A w 1ii(2Ad/cO. Treating the dust 
particle as a spherical capacitor, its potential and charge 
are related by <£ d = q/C with C = 27re <i(l + d/(2X D )). 

Evaluating ^ requires knowledge about ion density 
and velocities in the sheath. These quantities are di- 
rectly accessible by PIC and are shown in the lower panel 
of Fig. [9] Alternatively, the sheath parameters can be es- 
timated from Child's law^. The large dc-bias applied to 
the AE pixel locally dominates the characteristics of the 
sheath and masks the rf-character of the global discharge 
there. Low pressure should permit the neglect of ion col- 
lisions in the sheath and due to the additional bias the 
sheath potential is much larger than the electron tem- 
perature. Yet, within Child's law the sheath width is 
less than half the width of the PIC solution, with sim- 
ilar density and velocity profiles on this reduced scale. 
We attribute this discrepancy to the failure of Child's 
law to account for the local disturbance of the plasma 
conditions by the pixel biasing. Using the PIC data, 
we may calculate the two contributions to the ion drag 
force for different particle charges (see upper right panel 
of Fig. [9]). The bias switching instantaneously enlarges 
the sheath and thus the initial particle position is al- 
ready deep inside the sheath. Irrespective of the particle 
charge, the collection radius b c in this region essentially 
coincides with the particle radius. Only extremely high 
charges (q ~ — 10 5 eo) can compensate the high ion ve- 
locity resulting in a noticeable second term under the 
square root in Q. The proportionality of the Coulomb 
scattering radius b s to the particle charge often results in 
a dominance of b s over b c . In the relevant region, how- 
ever, the high values of V[ result in a small prefactor and 
b c and b s become comparable for q ~ —5 x 10 4 eo- Sum- 
ming up both contributions, the resulting ion drag force 
for q = — 10 5 eo is given in the upper left panel of Fig. [9] 
Due to the quadratic dependence of Fi on b c and 6 S , the 
curves for smaller q are marginal on this scale. In order 
to explain the remaining difference to the measured total 
force, the only explanation seems to be a positive recharg- 
ing of the particle. The resulting electric field forces for 
selected positive charges are given in left panel of Fig. [9] 
An explanation of the observed behavior requires a grad- 
ual charging to q ~ 2.5 x 10 4 eo during the fall. Cru- 
cial for this recharging procedure are: (i) absence of any 
electrons in the vicinity of the particle during the whole 
rf-cycle, which is guaranteed by the large negative pixel 
biasing; (ii) instantaneous sheath widening and location 
of the particle deep inside the widened sheath. Then 
the accelerated ions may acquire high enough energies to 
reach the particle surface, even for a positively charged 
particle. 
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IV. CONCLUSIONS 

In this work we studied the behavior of a micropar- 
ticle in the plasma sheath upon abrupt changes of the 
local plasma conditions. Thereto, we combine experi- 
mental studies with theoretical modeling and numerical 
simulations. Local perturbations of the plasma sheath 
in front of biased segments of the 'adaptive electrode' 
are monitored by the behavior of microparticles. Hereby, 
special emphasis is put on the particle dynamics upon 
bias switching and in the afterglow of the discharge. Re- 
laxation into a new equilibrium position is a markedly 
slower process than the adaption of the particle charge 
to the local electron and ion densities. Whether a par- 
ticle remains hovering above a biased segment or drops 
down depends on the charge equilibrium imposed by the 
altered plasma conditions in its vicinity. Slight changes 
of the plasma conditions initiate a damped oscillation of 
the particle into a new equilibrium position. Drastically 
increasing the negative bias of a pixel, locally the sheath 
may be deprived completely of electrons. This out-of- 
equilibrium situation may disturb the particle's charge 
balance to the point of a positive recharging. The re- 
sulting downward force might explain the observed addi- 
tional acceleration. Gravity and ion drag force alone are 
too weak to account for this behavior. 
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FIG. 10: (Color online) Left panel: Construction of image 
charges to calculate the force on a charged particle (black) 
between two parallel, grounded electrodes. Starting with the 
AE, an image charge of opposite sign (light grey) is con- 
structed at z = — zo. To describe the effect of the PE, image 
charges for both the original and the first image charge have 
to be constructed at z = 2d ± zq. These steps have to be 
iterated and generate an infinite series of image charges by 
repeatedly reflecting the newly created image charges at the 
AE and PE. The force acting on the original charge is then cal- 
culated from Coulomb's law accumulating contributions from 
all charges. Right panel: Force on a particle between two 
parallel grounded plates as a function of distance from the 
electrode. 



of magnitude smaller than the gravitational force acting 
on the considered particle. Thus, F s is ne gligib le except 
for the vicinity of the electrodes at which (A2) diverges. 



For the above parameters gravity and F s become compa- 
rable for zq ~ 0.05 mm. 



Appendix B: Details on PIC simulation 



Appendix A: Particle self-force 

Summing up the contributions of all image charges in 
Fig. [lOj we get 



F R = 



47Te 



1 

Az 1 



£ 

3=0 



(-iy 



(Al) 



with rj = (-l)Li/2j ( 2j L(Lj'/4J + 1) + (-lp' +1 ^o), where 
the floor function |_.J returns the largest integer that is 
not greater th an i ts argument, e.g., |_1 -75J = 1. The 
infinite sum in ( Al ) may be evaluated and with x — zo/L 
and the trigamma function ipi(x) = Y^jLo V( x + j) 2 we 
get 



• cot(7nr)) — ^ipi(x) 



(A2) 



In the right-hand panel of Fig. 10 we show F s as a func- 
tion of distance from the AE. For the rel evan t parameters 
(L = 10 cm, q ~ 10 4 eo) the prefactor in (A2) is six orders 



All presented simulation results have been calculated 
using an electrostatic 2d(r, z)3v PIC code. Since simu- 
lating the whole reactor vessel is infeasible, we restricted 
the simulation volume to the relevant region above the 
AE (L z = 5 cm, r max = 6 cm). By suitable choice 
of the substitute discharge parameters we ensured bulk 
plasma conditions in accordance with the experimental 
ones, $ p = 21.8 V, T e = 1.1 eV, n e = 2.5 x 10 9 cm- 3 . 
While the obtained plasma potential and densities devi- 
ate from their target values less than 1%, they are slightly 
larger for the electron temperature. 

Adopting the general code structure and null- collision 
met ho cP^^ from the xpdp2 code suit e,^^ our code has 
been tailored to match the specific requirements: rota- 
tional symmetry, inclusion of the AE, individual parti- 
cle weighting, OpenMP-parallelization). The neglect of 
Coulomb collisions between charged species can be justi- 
fied by the low degree of ionization (n e /n\ r ~ 10 -5 ) for 
the considered experimental conditions. Thermal motion 
of the fixed background gas is taken fully into account for 
ion neutral collisions while being neglected for electron 
neutral collisions due to the large mass ratio and thus 
faint momentum transfer. The solution of the Poisson 
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equation is performed by means of Fourier accelerated 
cyclic reduction^ as implemented, e.g., in FISHPACK^. 
We ensured numerical stability of our results for the con- 
sidered grid sizes (n z x n r = 1000 x 1200) by using double 
precision variables for the potential and charge density. 
The constant grid spacing of A r = A z = 0.05 mm is 
smaller than half the Debye length. By choosing a time 
step of At = 2.83 x 10 -11 s we ensure uj p At < 1/5, where 
u p is the plasma frequency. Respecting these constraints, 
the PIC method provides reliable results. 27 Modeling a 
vertical cut through the reactor vessel, symmetry implies 
a vanishing electric field at the inner boundary, r = 0, of 
the simulation volume. At the outer boundary, r = r max , 
we assume vanishing electric field boundary conditions 
since the lateral dimension of the simulation volume cov- 
ers only a fraction of the powered electrode. In order to 
model the AE, each pixel has been associated to an indi- 
vidual external circuit, and the charge balance is calcu- 
lated in each iteration. A linear weighting scheme for par- 
ticle deposition to the grid and interpolation of the elec- 
tric field ensures conservation of momentum and absence 
of self-forces. The r-dependence of the volume elements 
in cylindrical coordinates leads to a drastically varying 



number of superparticles in each cell when describing a 
constant density. Alternatively, additional weighting fac- 
tors for each superparticle may compensate for this^. 
Thereby, the statistics can be markedly improved since 
the initial number of superparticles in all cells is bal- 
anced. Throughout this work we used a mean number of 
superparticles per cell of (N c ) ~ 80. Note that these indi- 
vidual weighting factors have to remain constant during 
the whole simulation in order to ensure conservation of 
energy and momentum and absence of self-forces. As a 
consequence of the experimental setup, the main particle 
motion is in z-direction. Radially, the particles remain 
within nearby cells of their initial position. Thus, we are 
in the lucky position that the statistical advantage due 
to the individual weighting decays only slowly with time. 

In the experiment the AE-biasing does not influence 
the global discharge parameters but only locally disturbs 
the plasma sheath. This behavior reflects the small ratio 
of biased electrode to wall surface. The restricted size of 
the simulation volume implies a more pronounced impact 
of the AE-pixels on the bulk plasma due to the larger 
surface ratio. Nevertheless, the numerical bulk quantities 
with and without AE-biasing differ only by around 1%. 
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